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Abstract 

The existence of a light dark force mediator has been tested with the KLOE 
detector at DA$NE. This particle, called U, is searched for using the decay chain 
(j) ^ r]U, T] ^ vr+vr^vr'', U e~^e~. No evidence is found in 1.5 fb~^ of data. 
The resulting exclusion plot covers the mass range 5 < Mu < 470 MeV, setting 
an upper limit on the ratio between the U boson coupling constant and the fine 
structure constant, a' /a, of < 2 x IQ-^ at 90% C.L. for 50 < Mu < 420 MeV. 

Key words: e^e~ collisions, dark forces, gauge vector boson 
PACS: 14.70.Pw 



1 Introduction 



In recent years, several unexpected astrophysical observations have failed to 
find a common interpretation in terms of standard astrophysical or particle 
physics sources. A non-exhaustive list of these observations includes the 511 
keV gamma-ray signal from the galactic center observed by the INTEGRAL 
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satellite [T], the excess in the cosmic ray positrons reported by PAMELA [2], 
the total electron and positron flux measured by ATIC |3], Fermi [1], and 
HESS j5|6], the annual modulation of the DAMA/LIBRA signal [HH] and the 
low energy spectrum of nuclear recoil candidate events observed by CoGeNT 



Although there are alternative explanations for some of these anomalies, they 
could be all explained with the existence of a dark matter weakly interacting 
massive particle, WIMP, belonging to a secluded gauge sector under which the 
Standard Model (SM) particles are uncharged [TU|llfl2fl3fl4fl5fl(jfl7fl8fl9] . 
An abelian gauge field, the U boson with mass near the GeV scale, couples the 
secluded sector to the SM through its kinetic mixing with the SM hyper-charge 
gauge field. The kinetic mixing parameter, e, is expected to be of the order 
10~^-10~^ prpP] . so that observable effects can be induced in C(GeV)-energy 
e"^e~ colliders [2(J|21|22|23|24] and fixed target experiments [25|26|27|I25] . The 
possible existence of a new light boson gauging a new symmetry with a small 
coupling was in fact already introduced on general grounds in [29], and re- 
discussed in models postulating also the existence of light spin or 1/2 dark 
matter particles [SOfST] . This boson can have both vector and axial- vector cou- 
plings to quark and leptons, however axial couplings are strongly constrained 
by data, leaving room to vector couplings only. 

The U boson can be produced at e~^e~ colliders via different processes: e~^e^ — )■ 
U'y, e~^e~ — )■ Uh' (/I'-strahlung), where h' is a higgs-like particle responsible 
for the breaking of the hidden symmetry, and V P'-f decays, where V and 
P are vector and pseudoscalar mesons, respectively. In this work, we study 
the process (p ^ rjU, using a sample of (p mesons produced resonantly at the 
DA$NE collider. The U boson can be observed by its decay into a lepton pair, 
while the 1] can be tagged by one of its main decays. An irreducible background 
due to the Dalitz decay of the (p meson, cf) — )■ rii^i~, is present. This decay 
has been studied by the SND and CMD-2 experiments, which measured a 
branching fraction of BR{(j) r^e+e") = (1.19 ± 0.19 ± 0.07) x 10"^ and 
BR{(p r/e+e") = (1.14 ± 0.10 ± 0.06) x 10-^ respectively [32)133] . This 
corresponds to a cross section of cr(0 — )■ rji^i^) ~ 0.7 nb, with a di-lepton 
mass range < 470 MeV. For the signal, the expected cross section is 
expressed by [22] : 



where F^^{m\j) is the cpr]'^* transition form factor evaluated at the U mass 
while the following term represents the ratio of the kinematic functions of 



^ X^/^ {ml, ml) 
A3/2(m^,m2,0) 
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A(m^, 771,2, mg) ~ [1 + 7T7|/(777^ — m. 



1)]^ — 4777^7773/(777^ — 7773)^. 
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cr(0 — )■ ?7 f/) ~ 40 fb is obtained. Despite the small ratio between the overall 
cross section of (p ^ rjU and (p ^ rj , their different di-lepton invariant 
mass distributions allow to test the e parameter down to 10~^ with the KLOE 
data set. 

The best U decay channel to search for the (p ^ r]U process at KLOE is in 
e~^e~ , since a wider range of U boson masses can be tested and are easily 
identified using a time-of-flight (ToF) technique. The r] can be tagged by the 
three-pion or two-photon final state, which represent ~ 85% of the total decay 
rate. We have used the rj — ?► tt+tt^tt^ decay channel, which provides a clean 
final state with four charged particles and two photons. 



2 The KLOE detector 

The KLOE experiment operated from 2000 to 2006 at DA$NE , the Frascati 
0-factory. DA$NE is an e~^e~ collider running at a center-of-mass energy of 
~ 1020 MeV, the mass of the meson. Equal energy positron and electron 
beams collide at an angle of 7r-25 mrad, producing mesons nearly at rest. The 
detector consists of a large cylindrical Drift Chamber (DC), surrounded by a 
lead-scintillating fiber electromagnetic calorimeter (EMC). A superconducting 
coil around the EMC provides a 0.52 T field. The beam pipe at the interaction 
region is spherical in shape with 10 cm radius, it is made of a Beryllium- 
Aluminum alloy of 0.5 mm thickness. Low beta quadrupoles are located at 
about ±50 cm distance from the interaction region. The drift chamber pi], 
4 m in diameter and 3.3 m long, has 12,582 all-stereo tungsten sense wires and 
37,746 aluminum field wires. The chamber shell is made of carbon fiber-epoxy 
composite with an internal wall of ~ 1 mm thickness, the gas used is a 90% 
helium, 10% isobutane mixture. The spatial resolutions are a^y ~ 150 fim 
and (T2 ~ 2 mm. The momentum resolution is o'{pj_)/p± ~ 0.4%. Vertexes 
are reconstructed with a spatial resolution of ~ 3 mm. The calorimeter |35j is 
divided into a barrel and two endcaps, for a total of 88 modules, and covers 98% 
of the solid angle. The modules are read out at both ends by photomultipliers, 
both in amplitude and time. The readout granularity is ~ (4.4 x 4.4) cm^, for 
a total of 2440 cells arranged in five layers. The energy deposits are obtained 
from the signal amplitude while the arrival times and the particles positions are 
obtained from the time differences. Cells close in time and space are grouped 
into energy clusters. The cluster energy E is the sum of the cell energies. The 
cluster time T and position R are energy-weighted averages. Energy and time 
resolutions are as/E = h.l%/^E (GeV) and at = 57 ps/^E (GeV) © 100 ps, 
respectively. The trigger [36] uses both calorimeter and chamber information. 
In this analysis the events are selected by the calorimeter trigger, requiring 
two energy deposits with E > 50 MeV for the barrel and E > 150 MeV for 
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the endcaps. A cosmic veto rejects events with at least two energy deposits 
above 30 MeV in the outermost calorimeter layer. Data are then analyzed by 
an event classification filter |37j, which selects and streams various categories 
of events in different output files. 



3 Event selection 

The analysis of the decay chain — )■ 77 [/, — t- tt+tt^tt'^, U — )■ e+e^, has 
been performed on a data sample of 1.5 fb^'^, corresponding approximately 
to 5 X 10^ produced mesons. The Monte Carlo (MC) simulation of the 
irreducible background — )• rje^e', t] — )■ '7r"'"7r~7r° has been produced with 
dr{(j) — )■ 7] e~) / dm^e weighted according to Vector Meson Dominance model 
|38] . using the form factor parametrization from the SND experiment [32] . 
The MC simulation for the (p ^ rjU decay has been developed according to 
[22] . with a fiat distribution in M^e- All MC productions, including all other 
(j) decays, take into account changes in DA$NE operation and background 
conditions on a run-by-run basis. Data-MC corrections for cluster energies 
and tracking efficiency, evaluated with radiative Bhabha events and — ?■ pvr 
samples, respectively, have been applied. 

As a first step of the analysis, a preselection is performed requiring: 

(1) two positive and two negative tracks with point of closest approach to 
the beam line inside a cylinder around the interaction point (IP), with 
transverse radius Rpy = 4 cm and length Zpy = 20 cm; 

(2) two photon candidates i.e. two energy clusters with E > 7 MeV not 
associated to any track, in an angular acceptance | cos 6*^1 < 0.92 and in 
the expected time window for a photon (|T^ — i?^/c| < MIN(5(Jt, 2 ns)); 

(3) best 7r"'"7r~77 match to the 7] mass in the pion hypothesis to assign tt^ 
tracks[ll the other two tracks are then assigned to e^; 

(4) loose cuts of about ± 4 a's on 77 and vr" invariant masses (495 < M^+^-^^ < 
600 MeV, 70 < M^^ < 200 MeV). 

After this selection, a clear peak corresponding to — )■ 77 e"'"e~ events is ob- 
served in the distribution of the recoil mass to the e~^e~ pair, M^ecoi\{ee), as 
shown in Fig. [TJ The second peak at ~ 590 MeV is due to Ks — ?• n'^ir' events 
with a wrong mass assignment. Events in the 535 < Mj.ecoii(ee) < 560 MeV 
window are retained for further analysis. 

A residual background contamination, due to — )■ 777 events with photon 

^ The invariant mass of 7r"'"7r~77 for each positive/negative track pair, Mtestj is 
evaluated in the hypothesis that the two tracks belong to charged pions. The track 
pair with the smaller |Mtest — Af?j| value is assigned to and tt~ . 
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Fig. 1. Recoiling mass against the e^e~ pair for the data sample after preselection 
cuts. The (p r]e~^e~ signal is clearly visible in the peak corresponding to rj mass. 
The second peak at ~ 590 MeV is due to — )• KsKl, Kg — > vr^vr" events with 
wrong mass assignment. 



conversion on beam pipe (BP) or drift chamber walls (DCW), is rejected by 
tracing back the tracks of the two e"^, e~ candidates and reconstructing their 
invariant mass {M^e) and distance (-Dee) at the BP/DCW surfaces. As both 
quantities are small in case of photon conversions, — )■ 777 background is 
removed by rejecting events with: Mee(-BP) < 10 MeV and Dee{BP) < 2 
cm, M^^{DCW) < 80 MeV and D^^iDCW) < 10 cm. A further relevant 
background, originated by — KK decays surviving analysis cuts, has more 
than two charged pions in the final state and is suppressed using time-of-flight 
(ToF) to the calorimeter. When an energy cluster is connected to a track, the 
arrival time to the calorimeter is evaluated both using the calorimeter timing 
(^cluster) and the track trajectory (Tt^ack = -f/track/z^c). The AT = 

-^track -^cluster 

variable is then evaluated in both electron (ATe) and pion (AT^) hypotheses. 
Events with an e"*", e~ candidate outside a Sa's window on the ATe variables 
are rejected. In Fig. El the M^e distribution evaluated at different steps of 
the analysis is shown. The peaks at ~ 30 MeV and ~ 80 MeV are due to 
photon conversions on BP and DCW, respectively. The ToF cut reduces the 
tail at high Mee values while the conversion cut removes events in the low 
invariant mass region. The analysis efficiency, defined as the ratio between the 
number of events surviving analysis cuts and that of all generated events, is 
shown in Fig. |3] as a function of Mee, ranging between 10% and 20%. The 
main contribution to the loss of efficiency is due to preselection cuts, being 
epresel = (24.73 ± 0.04)%. 

In Fig. m the comparison between data and Monte Carlo events for Mee and 
cosip* distributions is shown. Here ip* is the angle between the r] and the e"*" 
in the e'^e~ rest frame. About 14,000 — )■ rie~^e~, rj — )■ -7r"'"7r~7r° candidates 
are present in the analyzed data set, with a negligible residual background 
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Fig. 2. Alee distribution for data after different analysis cuts. 
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Fig. 3. Analysis efficiency as a function of e~^e~ invariant mass. 



contamination. 



4 Upper limit evaluation 



As an accurate description of the background is crucial for tlie searcli of the 
U boson, its shape is extracted directly from our data. A fit is performed 
to the Mee distribution, after applying a bin-by-bin subtraction of the — )■ 
7^7 background and efficiency correction. The parametrization of the fitting 
function has been taken from Ref. j38] : 
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Fig. 4. Invariant mass of the e~^e~ pair (top) and cosV'* distribution (bottom) for 
(j) r]e^e~, rj — )> 'k~^t:~'k^ events. Dots are data, the black solid line is the sum 
of all MC expectations while signal and residual background contamination from 
(/) — 7> 777 are shown in colors. 



(ir(0 — )■ e'^e 



Am? 



2m' 



\ A3/2( 



(2) 



with q = Mee and the transition form factor described by 



Free parameters of the fit are A and an overall normalization factor. A good 
description of the Mee shape is obtained except at the high end of the spec- 
trum (see Fig. [5l), where a residual background contamination from multi-pion 
events is still present. 

As mentioned in Sec. [3], the (p ^ rjU MC signal has been produced according 
to Ref. [22], with a flat distribution of the U boson invariant mass, M[/. This 
sample has been used to evaluate the resolution on the e~^e~ invariant mass 
as a function of Mu, applying a Gaussian fit to the Mee — distributions. 
Results are reported in Fig. [61 The resolution is ~ 2 MeV for Mu < 350 MeV 
and then improves to 1 MeV for higher values. The upper limit on — )■ 77 [/ 
signal as a function of Mu is then obtained in the following way: 
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Fig. 5. Fit to the corrected Mge spectrum for the Dalitz decays (/> — ?■ rje~^e~ . 

(a) MC events are divided in sub-samples of 1 MeV width in the range 5 < 
Mu < 470 MeV; 

(b) for each Mjj sub-sample, the average value of the (j) — )■ rje^e~ background, 
b{Mee), is obtained by fitting the reconstructed Mee spectrum with 5 MeV 
binning, removing five bins centered at Mjj; 

(c) for each fit, the maximum variation of b{Mee) events, A6(Mee), is obtained 
changing by ±1 a the fit parameters; 

(d) for each Mu value, the signal hypothesis is tested comparing observed data, 
b{Mee) and MC signal in the five reconstructed bins excluded in (b). The 
exclusion plot is obtained applying the CLs method |39] . A Gaussian spread 
of width A6(Mee) on the background distribution is applied while evaluating 
CLs. 

In Fig. [7]the exclusion plot at 90% C.L. on the number of events for the decay 
chain (p ^ rjU, r] ^ 7r~^'ir~'ir^, U — > e^e~, is shown. Using Eq. ([T]) and taking 
into account the analysis efficiency this result is then reported in terms of the 
parameter a' /a = e^, where a' is the coupling of the U boson to electrons 
and a is the fine structure constant. The opening of the U — > fj^^fi~ threshold, 
in the hypothesis that the U boson decays only to lepton pairs and assuming 
equal coupling to e~^e~ and fi'^fJ^', has been included. In Fig. [8] the smoothed 
exclusion plot at 90% C.L. on a' /a is compared with existing limits from the 
muon anomalous magnetic moment [10] and from recent measurements of 
the MAMI/Al [H] and APEX |12] experiments. The gray line is where the U 
boson parameters should lay to account for the observed discrepancy between 
measured and calculated values. Our result greatly improves existing limits 
in a wide mass range, resulting in an upper limit on the a' /a parameter of 
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Fig. 6. Resolution on Mf>e as a function of the U boson invariant mass for (p ^ rjU 
MC events. 
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Fig. 7. Upper limit at 90% C.L. on the number of events for the decay chain (p rjU, 
rj i:~^TT~-K^ , U — >■ e^e~. 

< 2 X 10"^ @ 90% C.L. for 50 < Mc/ < 420 MeV, thus covering part of the 
expected e range (see Sec. [1]). We exclude that the existing discrepancy is 
due to a ?7 boson with mass ranging between 90 and 450 MeV. 
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